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OBJECTIVE 

Develop  and  validate  a  technique  for  passively  inferring  the 
tropospheric  refractivity  structure  from  observations  of  low- 
angle,  satellite- to-ground  rf  transmissions. 

RESULTS 

1.  Eleven  measurements  of  dual-frequency  radio  transmissions 
were  made  during  the  summer  of  1978. 

2.  Application  of  the  technique  correctly  infers  the  ducting 
structure  in  seven  of  the  lower  frequency  observations  and  five 
of  the  high-frequency  observations. 

3.  Inferred  profile  structure  can  differ  significantly  from 
the  observed  refractive  structure. 

4.  The  major  limitation  to  the  applicability  of  this  techni¬ 
que  is  an  uncorrectable  rms  range  error  of  approximately  4  km. 

5.  Moderate  success  was  achieved  in  inferring  the  refractive 
structure  from  rf  measurements. 

RECOMMENDATION 

Terminate  this  project,  since  operationally  significant  data 
cannot  be  extracted  reliably. 
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INTRODUCTION 


Modern  naval  microwave  communication  and  sensor  system  per¬ 
formance  can  be  greatly  affected  by  refractive  conditions-  For 
example,  a  typical  surface  search  radar  system,  under  standard 

atmospheric  conditions,  can  normally  detect  targets  up  to  30  km 
distant.  Under  atmospheric  conditions  in  which  a  surface-based 
duct  is  present,  however,  it  is  not  uncommon  for  this  radar 
system  to  detect  the  same  targets  at  ranges  in  excess  of  300  km. 

Statistics  on  refractive  conditions  found  in  a  maritime 

environment  are  useful  as  an  aid  in  determining  the  scope  of 
environmental  effects  on  microwave  system  performance.  A  5-year 
data  base  of  upper  air  soundings  from  921  radiosonde  stations 
distributed  throughout  the  world  was  extensively  analyzed  by 
GTE/Sylvania  (ref  1).  The  results  of  this  study  were  examined  by 
NOSC  for  radiosonde  stations  situated  along  ocean  coastal 
regions,  island  stations,  and  meteorological  station  ships  (cate¬ 
gory  "C"  stations).  The  following  conclusions  were  drawn:  (1) 

Surface-based  ducts  occur  worldwide  6  percent  of  the  time  on  the 

average.  (2)  The  average  thickness  of  the  surface-based  duct  is 
less  than  100  m.  (3)  Elevated  ducts  occur  13  percent  of  the  time 
on  the  average.  (4)  The  average  height  of  the  elevated  duct  is 
less  than  1500  m.  These  percentages  seem  low,  but  they  represent 
an  average  worldwide  occurrence.  In  some  areas  of  the  world,  the 
percentages  are  significantly  higher.  In  the  San  Diego  locale, 
for  example,  a  surface-based  duct  occurs  23  percent  of  the  time 
and  an  elevated  duct  occurs  48  percent  of  the  time.  These 
percentages  relate  directly  to  microwave  system  performance  and 
can  be  used  to  indicate,  by  geographical  regions,  locales  in 
which  propagation  anomalies  rr  «y  significantly  affect  operations. 
While  such  statistics  may  be  useful  for  long-range  planning,  it 
is  the  current  refractive  couditions  that  must  be  understood  for 
immed ia te  operations. 

Assessment  of  the  refractive  effects  on  microwave  systems 
must  be  considered  a  twofold  problem.  First,  a  reasonably 
accurate  representation  of  the  vertical  refractivity  structure  is 
required.  Second,  a  computational  and  display  system  is  needed 
that  relates  the  refractive  environment  to  those  effects  by  means 
of  EM  propagation  models.  The  Integrated  Refractive  Effects 
Prediction  System  (IREPS)  (ref  2)  provides  the  hardware  and 
software  necessary  to  determine  the  atmospheric  refractive 
effects  on  microwave  systems,  IREPS  is  currently  installed 
aboard  13  of  the  Navy's  major  combatants  and  its  products  are 
valued  in  operational  decisionmaking.  This  system  essentially 
satisfies  the  second  portion  of  the  assessment  problem. 


1.  GTE/Sylvania  Inc  Final  Report,  Radiosonde  Data  Analysis  Five 
Years  Data,  by  LD  Ortenburger,  NSA  Contract  MDA  904-76-C-0233  ,  29 
July  1977. 

2.  NOSC  TD  238,  Integrated  Refractive  Effects  Prediction  System 
(IREPS),  Interim  User's  Manual,  by  HV  Hitney  and  RA  Paulus,  March 
1979  . 
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Two  methods  are  used  by  th^  operational  forces  to  obtain  the 
refractivity  structure;  a  balloonborne  radiosonde  and  an  air¬ 
borne  microwave  ref rac tome ter .  Each  method  suffers  drawbacks.  In 
particular,  the  radiosonde  radiates  approximately  1  W  at  403  MHz, 
which  renders  this  method  undesirable  during  periods  of  electro¬ 
magnetic  control  ( EMCON ) .  Because  it  is  an  airborne  system,  the 
ref ractome ter  is  an  expensive  measurement  and  is  available  only 
on  a  select  few  aircraft* 


Although  each  method  can  measure  the  refractivity  structure, 
a  technique  entitled  Passive  Refractive  Index  by  Satellite 
Monitoring  (PRISM)  was  proposed  as  a  potentially  useful  third 
alternative.  By  receiving  rf  transmissions  from  a  satellite  and 
inferring  the  refractivity  structure  from  the  location  of  the 
interference  nulls  (created  by  the  path-length  difference  between 
the  direct  and  sea-re£lected  rays) ,  PRISM  provides  several  advan¬ 
tages.  First  and  foremost,  the  technique  is  completely  passive 
since  only  the  satellite  radiates.  Second,  since  the  location  of 
the  interference  null  is  detectable  by  measuring  the  relative 
amplitude  of  the  signal,  the  receiver  hardware  is  relatively 
simple.  Third,  the  satellite  does  not  have  to  be  tracked  by 
ground- station  equipment  since  only  a  useable  signal  is  required. 
Lastly,  logistics  support  is  minimal  compared  to  that  of  radio¬ 
sonde  and  ref rac tome  ter  methods  . 


It  is  the  purpose  of  this  paper  to  present  the  background, 
procedures,  measurements,  and  results  of  the  PRISM  Project, 
PRISM  was  initiated  in  October  1977  under  the  sponsorship  of  the 
Naval  Sea  Systems  Command,  with  the  objective  to  develop  and 
validate  a  technique  proposed  by  HV  Hitney  (ref  3)  to  infer  the 
tropospheric  refractivity  structure  passively  from  rf  measure¬ 
ments. 


3.  NOSC  Technical  Note  NELC  TN  3108,  A  Proposed  Method 
Refractive-Index  Profiles  by  Satellite  Monitoring,  by  HV 
21  January  1976  .  NOSC  TN  s  are  informal  documents 
chiefly  for  internal  use. 


to  Infer 
Hitney , 
intended 
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BACKGROUND 


Many  experiments  have  been  performed  in  the  general  field  of 


rf  inversion  (ref  4 
ful  in  determining 


Occultation  experiments  have  been  success- 
the  refractivity  structure  for  altitudes 
generally  in  excess  of  10  km.  Less  successful  have  been  attempts 
to  infer  the  structure  from  measurements  of  elevation  angle  and 
Doppler  shift.  These  methods  generally  require  either  extremely 
high  accuracy  of  the  ground-station  tracking  equipment  or  high 
frequency  stability.  In  1977,  when  this  study  began,  no  inver¬ 
sion  technique  was  available  to  infer  the  low-level  (less  than  4 
km  altitude)  refractivity  structure  passively  with  significant 
reliability  and  precision. 

Examination  of  radio  data  measured  on  an  aircraft- to-ground 
path  (ref  5)  indicated  that  the  location  of  the  interference 
nulls  could  be  used  first  to  infer  the  elevation  angles  to  the 
required  accuracy,  then  to  infer  the  refractivity  profile.  In 


figure  1, 
elevation 

the  path  length  difference  6  is  related  to 
angle  a  (by  ref  6)  as  follows: 

the 

true 

9  = 

/(a/3)2  +  2h/3ae  -  a/3 

a) 

6  = 

2  (a  +  6)2  /h2  +  a  (a  +  h)  92 
v  e  e 

(2) 

Equations  (1) 
are  parallel, 
the  aircraft, 
profiles  were 

and  (2)  are  valid  if  the  direct  and  reflected  rays 
Considering  the  necessarily  limited  altitude  of 
this  condition  was  not  strictly  met.  However, 

inferred  that  closely  approximated  the  observed, 

warranting  the 
ground  path  on 


extension  of  this  technique  to 
which  the  rays  would  be  parallel. 


a  satellite- to- 


A  ground-station  receiver  was  designed  and  a  suitable  satel¬ 
lite  was  chosen  in  the  time  frame  from  October  1977  to  June  1978. 
During  July  and  August  of  1978,  observations  of  dual  rf  transmis¬ 
sions  and  the  local  refractivity  structure  were  made  (ref  7). 


4.  NASA  TM  X- 62,  150,  Mathematics  of  Profile  Inversion,  edited  by 
L  Colin,  August  1972. 

5.  NOSC  Technical  Note  NELC  TN  2874,  Propagation  Measurements  at 
3  GHz  in  the  San  Diego  Offshore  Area,  by  KD  Anderson  and  DR 
Jensen,  20  January  1975. 

6.  NRL  Report  7098,  Machine  Plotting  of  Radio/Radar  Vertical 
Coverage  Diagrams,  by  LV  Blake,  Naval  Research  Laboratory,  25 
June  1970  . 

7.  NOSC  TN  639,  PRISM  —  Passive  Refractive  Index  by  Satellite 

Monitoring:  Preliminary  Measurements  and  Results,  by  KD  Anderson 

and  HV  Hitney,  March  1979. 
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Application  of  this  inference  technique  to  the  observed  data 
illustrated  a  major  problem  in  the  determination  of  6:  the  path- 
length  difference  could  not  be  extracted  reliably  from  the  fre¬ 
quency  diversity  information  available.  A  procedure  to  recover  6 
from  a  simplified  atmospheric  model  (ref  8)  was  investigated. 
This  procedure  started  with  higher  elevation  angles  (70  milli- 
radians) ,  where  the  refractivity  structure  has  negligible  effect, 
and  iterated  6  for  the  lower  angles.  But  since  the  combination 
of  measurement  inaccuracies  and  low  angle  refraction  effects 
dominated,  this  procedure  was  found  to  be  untenable. 

Since  the  extraction  of  <$  proved  unworkable,  the  inference 
technique  was  modified.  Briefly,  it  was  found  that  possible  null 
locations  are  easily  determined  from  the  signal  received.  Some 
of  these  locations  must  be  true  nulls.  Comparing  all  possible 
null  locations  found  to  a  library  of  null  locations  (computed  by 
ray-trace  techniques  for  a  wide  selection  of  assumed  refractivity 
profiles)  could  yield  a  respectable  inferred  profile. 

The  following  sections  of  this  paper  will  examine  the 
methodology  of  the  modified  technique,  the  measurements,  and 
results. 


EXPERIMENTAL  PROCEDURES 


METHODOLOGY 

Inference  of  the  tropospheric  refractivity  ducting  environ¬ 
ment  is  based  upon  a  statistical  analysis  of  satell ite- to-ground 
rf  transmissions.  When  the  satellite  is  near  the  horizon  (fig 
1),  the  received  signal  is  the  sum  of  two  rays.  One  ray  propa¬ 
gates  only  through  the  atmosphere,  the  direct  path,  while  the 
second  ray  is  reflected  from  the  surface  of  the  ocean.  Since 
there  is  a  phase  shift  of  n  radians  upon  reflection,  the  received 
signal  is  a  maximum  when  the  path  phase  difference  between  the 
direct  and  reflected  ray  is  an  integral  multiple  of  it.  Conver¬ 
sely,  when  the  path  phase  difference  is  a  multiple  of  2 tt  (i.e. 
the  path  length  difference  is  a  multiple  of  one  wavelength) ,  the 
received  signal  is  a  minimum.  As  the  satellite  moves  through  its 
orbit,  the  path  length  difference  changes  over  many  wavelengths, 
resulting  in  a  received  signal  that  repeats  maxima  and  minima 
with  time  or,  correspondingly,  with  range.  Each  occurrence  of 
the  minimum,  or  null,  is  assigned  a  number  based  on  the  integral 
number  of  wavelengths  difference  between  the  direct  and  the  sea- 
reflected  paths.  Null  1  corresponds  to  a  path  length  difference 
of  one  wavelength,  null  2  to  a  difference  of  two  wavelengths,  and 
so  on. 


8.  A  Null  Finding  Algorithm  for  PRISM,  prepared  for  NOSC  by  PH 
Levine  et  al  under  Contract  N00123-78-C-0043 ,  23  July  1979. 
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Figure  2  portrays  the  form  of  the  received  signal  plotted  as 
amplitude  in  dB  versus  ground  range  from  the  receiver  to  the 
satellite.  Two  sets  of  curves  are  shown,  illustrating  the 
effects  of  two  different  refractivity  profiles.  The  solid-line 
plot  represents  the  signal  received  under  standard  atmospheric 
conditions  whereas  the  broken  line  shows  the  pattern  when  a  mod¬ 
erately  strong  surface-based  duct  is  present.  Both  cases  are 
computed  from  ray  trace  procedures  for  a  frequency  of  1239  MHz, 
satellite  height  of  1000  km  and  a  receiver  height  of  35  m.  In 
addition  to  the  received  signal  amplitude  curves,  the  two  curves 
proceeding  from  the  lower  left  to  the  upper  right  relate  the 
respective  path  length  differences  (null  numbers)  to  ground 
range.  The  first  null  observed  for  the  standard  profile  (solid¬ 
line  plot)  corresponds  to  a  path-length  difference  of  one 
wavelength.  It  occurs  at  a  range  of  3415  km.  Note  that  for  the 
surface-based  duct  environment  (broken-line  plot)  the  first 
observable  minimum  is  not  null  1,  but  rather  null  3,  which  is 
found  at  a  range  of  3405  km.  Note  also  that  as  the  path-length 
difference  increases  (at  increasing  elevation  angles)  the  loca¬ 
tion  of  each  null  tends  to  become  independent  of  the  refractivity 
structure . 

These  two  theoretical  plots  of  path-length  difference  as  a 
function  of  range  illustrate  the  methodology  of  the  inversion 
technique.  Simply,  the  form  of  the  path-length  difference  (pld) 
curve  is  extracted  from  the  observed  data  by  noting  the  range  at 
which  each  minimum  occurs,  without  regard  to  an  absolute  null 
number  assignment.  The  observed  pld  curve  is  then  compared  to  a 
library  of  pld  curves  generated  from  a  wide  variety  of  refrac¬ 
tivity  profiles.  The  inferred  profile  is  selected  as  the  refrac¬ 
tivity  structure  that  yields  the  minimum  standard  deviation  of 
the  difference  between  the  predicted  and  the  observed  ranges.  For 
a  given  library  profile  null  number,  the  range  at  which  the  null 
occurs  is  compared  to  every  possible  null  location  found  in  the 
observed  data.  If  the  minimum  difference  is  less  than  10  km,  the 
observed  location  is  noted.  Otherwise  the  library  null  number  is 
tagged  as  rejected.  The  minimum  range  difference  is  found  for 
each  null  of  the  library  profile.  Up  to  20  predicted  nulls  are 
examined.  The  standard  deviation  of  the  range  differences  is 
computed  for  the  library  profile.  After  these  procedures  have 
been  completed  for  every  profile,  the  "best  fit"  of  the  observed 
to  the  calculated  pld  curves  is  given  by  the  profile  that  has  the 
minimum  standard  deviation  of  the  range  differences  and  the 
fewest  number  of  rejected  library  nulls. 

Since  the  inference  technique  relies  upon  comparing  modeled 
to  observed  data,  the  prediction  accuracy  is  governed,  in  part, 
by  the  adequacy  (number  and  scope)  of  the  selected  library  pro¬ 
files.  Accordingly,  figure  3  is  a  plot  of  null  number  as 
ordinate  and  ground  range  as  abscissa  for  all  the  438  refractiv¬ 
ity  profiles  contained  in  the  library.  The  figure  illustrates 
the  pld  for  a  transmission  frequency  of  1239  MHz,  a  satellite 
height  of  1000  km  and  a  receiver  height  of  35  metres.  An  expan¬ 
sion  about  the  first  five  null  numbers  is  shown  as  figure  4.  On 
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this  figure  is  imposed  the  regional  influence  by  ducting  environ¬ 
ment  type.  That  is,  the  figure  illustrates  the  possible  ranges, 
as  a  function  of  pld,  that  can  occur  for  a  nonducting,  an 
elevated,  and  a  surface-based  ducting  refractivity  profile. 
There  are  regions  in  which  each  ducting  case  singularly  domin¬ 
ates,  and  there  are  regions  in  which  two  and  even  all  three 
ducting  cases  compete.  Although  similar,  the  various  pld  curves 
are  unique. 

Figures  5  and  6  show  the  received  signal  measured  on  20  July 
1978  at  1946  (GMT)  for  transmission  frequencies  of  1239  MHz  and 
2891  MHz,  satellite  height  of  1001  km,  and  receiver  height  of  35 
metres.  The  path  length  difference  curves  are  computed  from  ray 
trace  methods  by  using  the  refractivity  structure  as  measured  by 
a  radiosonde  taken  near  the  receiver  site  at  1945  (GMT)  on  the 
same  day.  As  expected,  the  interference  patterns  exhibit  maxima 
and  minima  in  the  presence  of  noise.  Noise  contributions  are 
particularly  noticeable  at  ranges  exceeding  3350  km. 

From  examination  of  the  data,  one  of  the  first  questions  to 
be  addressed  is  "What  constitutes  an  interference  null?"  For  the 
lower  frequency  data,  the  null  locations  for  ranges  less  than 
3350  km  are  obvious.  At  greater  ranges,  however,  are  the  minima 
true  nulls  or  noise?  Examination  of  the  higher  frequency  data 
indicates  fewer  obvious  nulls,  with  most  minima  uncertain.  This 
problem  was  examined  in  detail,  and  empirical  guidelines  were 
established  to  define  a  null.  These  rules  are  developed  as 
follows. 

A  range  window  is  centered  on  each  sample  point  of  the  data. 
If  the  amplitude  of  the  current  sample  is  the  minimum  value 
within  the  bounds  of  the  window,  the  range  at  which  the  current 
sample  occurs  is  tagged  as  a  possible  null  location.  The  width 
of  the  window  is  approximately  the  range  separation  between  adja¬ 
cent  high  order  nulls,  and  is  20  km  for  the  1239  MHz  data  and  10 
km  for  the  2891  MHz  data.  In  addition,  nulls  are  defined  to 
occur  at  ranges  less  than  the  maximum  range  at  which  the  received 
signal  is  greater  than  the  average  amplitude  of  all  samples.  The 
primary  motivation  for  this  second  requirement  is  explained  with 
reference  to  figures  5  and  6.  Note  that  the  signal  is  observed 
for  ranges  in  excess  of  3500  km  whereas  the  prediction  indicates 
that  the  first  observable  null  occurs  at  a  range  of  3409  km. 
Also  note  that  the  amplitude  of  the  signal  tends  to  decrease  at  a 
rate  of  0.1  dB  per  km  for  ranges  in  excess  of  3400  km.  Since  the 
measured  refractivity  structure  shows  a  strong  surface-based 
duct,  it  is  thought  that  the  unexpected  appearance  of  the  signal 
is  due  to  coupling  of  energy  into  the  duct  at  ranges  well  beyond 
the  horizon.  The  imposition  of  the  maximum  range  criteria 
removes  this  precursor  from  processing.  To  summarize  the  null 
definition,  a  maximum  range  at  which  nulls  occur  is  imposed  on 
the  basis  of  the  average  value  of  the  signal,  and  possible  null 
locations  are  selected  as  the  minimum  value  within  a  predefined 
range  window. 


The  following  sections  will  deal  with  the  development  and 

formulation  of  the  library  refractivity  profiles,  ray  tracing 
procedures,  and  the  hardware. 

REFRACTIVITY  PROFILE  MODEL 

The  refractivity  index,  n,  is  defined  as  the  ratio  of  the 
propagation  velocity  in  a  vacuum  to  the  velocity  in  the  medium. 

In  the  atmosphere,  the  value  of  n  differs  from  unity  by  about 

0.05  percent.  It  is  convenient  to  define  the  refractivity,  N,  as 
(n-1)  X  10  ,  which  has  values  in  the  range  0  to  500.  For  radio 

frequencies,  the  refractivity  is  related  to  atmospheric  tempera¬ 
ture,  pressure,  and  humidity  by  the  relationship 

N  =  77.6P/T  +  3.73  X  105e/T2,  (3) 

where  P  is  the  pressure  in  millibars,  T  is  the  air  temperature  in 
kelvins  and  e  is  the  saturated  water  vapor  pressure  in  millibars. 
In  some  applications  it  is  desirable  to  remove  the  effects  of  the 
earth's  curvature.  This  is  accomplished  by  the  modified  refrac¬ 
tivity,  M,  which  is  related  to  N  as  follows: 

M  =  N  +  0.1572,  (4) 


where  z  is  the  height  above  the  earth's  surface,  in  metres* 

Whenever  M  decreases  with  height,  a  so-called  trapping  layer 
is  formed.  An  EM  wave  entering  this  region  at  a  small  angle  is 
refracted  or  bent  at  a  curvature  greater  than  the  curvature  of 
the  earth's  surface.  The  height  interval  between  the  top  of  the 
trapping  layer  and  the  lowest  height  where  a  trapped  ray  propa¬ 
gates  is  called  the  ducting  region,  or  duct.  This  region  may  be 
found  by  simple  inspection  of  a  plot  of  M  as  abscissa  and  height 
as  ordinate  (fig  7).  A  straightedge  is  placed  parallel  to  the 
height  axis  through  the  M-unit  value  at  the  top  of  the  trapping 
layer,  and  the  ducting  region  is  the  height  interval  from  that 
point  to  the  first  lower  intersection  of  the  straightedge  with 
the  M  curve.  In  figure  7,  the  plot  of  M  vs  height  illustrates 
two  ducting  regions.  An  elevated  duct  exists  if  the  base  of  the 
duct  is  above  the  earth's  surface,  whereas  a  surface-based  duct 
exists  if  the  base  of  the  duct  is  at  the  earth's  surface. 


Empirically,  the  environment  model  was  divided  into  three 
regions  (fig  8).  Following  Bean  and  Dutton  (ref  9)  in  the  height 
interval  from  9  to  50  km  above  the  earth,  the  N  profile  varies 
exponentially  with  height  as  follows: 

N  =  105  exp[ -0.1424  ( Z-9 )  j ,  (5) 


where  Z  is  the  height  in  kilometres.  As  the  altitude  increases 
to  50  km,  N  becomes  0.31;  this  value  implies  that  the  medium  is  a 
near-vacuum.  For  heights  in  excess  of  50  km  the  medium  is 
assumed  to  be  a  vacuum  and  N  is  independent  of  height,  with  a 
value  of  zero.  In  the  region  from  the  surface  to  9  km  the  N 
curve  is  segmented,  varying  linearly  with  height.  On  the  basis 
of  many  observations  of  real-world  N  profiles  and  experience  in 
assessing  propagation  effects,  a  four-segment  construction  is 
used  for  the  linear  portion  of  the  library  N  profiles. 


RAY  TRACING  PROCEDURES 

The  location  of  the  minima  and  the  elevation  angle  in  the 
inversion  Library  are  derived  from  refractivity  profiles  previ¬ 
ously  described  and  ray  tracing  calculations  made  on  a  digital 
computer.  Since  these  computations  are  critical  to  the  inversion 
process  and  differ  from  procedures  developed  by  others,  this 
section  briefly  describes  the  formulation. 

A  spherical  coordinate  system  (p,  6,  p )  with  center  at  the 

center  of  the  earth  and  a  point  source  located  at  a  height  zi 

above  the  surface  of  the  earth  is  illustrated  in  figure  9.  Since 

the  refractivity  structure  is  assumed  to  be  spherically  symmet¬ 
ric,  all  functions  are  independent  of  p  and  only  the  meridian 
plane  is  considered.  Snell's  law  in  this  case  may  be  written  as 

pn  , ,  x 

—  cos  a  =  c,  (6) 

a 

where  a  is  the  earth's  radius,  p  is  the  vector  from  the  earth's 
center  to  a  point  along  the  ray,  a  is  the  angle  formed  between 

the  normal  of  p  and  the  ray,  and  c  is  a  constant  that  is  eval¬ 

uated  at  any  point  along  the  ray  where  n,  p,  and  a  are  known. 
From  simple  trigonometry,  eq  (6)  can  be  expressed  as 


pn 

—  sm  a 
a 


*/ 


( P n/a  ) 


T 

c  . 


(7) 


9.  Radio  Meteorology,  BR  Beam  and  EJ  Dutton,  p  64,  Dover 
Publications  Inc,  1968. 
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Combining  eq  (6)  and  (7), 


tan  ex  = 


r{  T72 - 2 

±M_pn 


/a)  -  c  =  dp 

c  pd  0  * 


Integrating  in  terms  of  the  central  angle  9, 


9  =  ±c  / 


dp 


p4pn/a)2  -  c2 


(8) 


(9) 


The  quantity  r,  which  is  the  distance  the  ray  travels  as  measured 
along  the  earth's  surface,  is  given  by  the  relationship 


r  =  0a . 


(10) 


Since  the  refractivity  structure  is  composed  of  three 
reg ions-1 inear ,  exponential,  and  constant  (unity)-the  solution  of 

eq  (9)  is  examined  in  each  region. 


LINEAR  REGION  RAY  TRACING 


In  the  height  interval  from  0  to  9  km  above  the  earth's 
surface,  it  is  assumed  that  the  distribution  of  n  can  be 
expressed  in  the  form  n(p)  =  bQ  +  for  an  interval  of  p.  The 

linear  variation  of  n  allows  a  closed  form  solution  to  eq  (9)  as 
follows: 


9 


-c 

/IT 


In 


2  m  +  p  t  + 
P 


(u  >  0) 


arc  tan 


2  m  +  p  t 
2  /v 


(u  <  0) 


(ID 


2/v 
c  - - 


(u  =  0) 


where 


/ 


t  =  2(b0bi)/a 
P  —  ( bg/ 3 ) 

v  =  l(bgp  +  b^)/a]2  -  c2 • 


IS 


For  each  linear  segment  of  ref ractivity ,  the  coefficients 
and  b-^  are  readily  determined.  A  computer  solution  to  this  equa¬ 
tion  is  easily  amenable  and  the  solution  is  rapidly  found. 

EXPONENTIAL  REGION  RAY  TRACING 

The  refractivity  in  the  height  interval  of  9  to  50  km  above 
the  earth's  surface  is  given  by  equation  (5)  as 

N  =  105  exp[-0.1424(Z-9)]  . 


For  any  known  pairs  of  a  and  n  at  the  surface,  the  range,  r,  in 
this  region  can  be  found  independently  from  all  other  regions 
since  this  curve  is  fixed.  If  the  surface  launch  angle,  clq,  of 
the  ray  and  the  surface  value  of  the  retractive  index,  Hq,  are 
given,  the  ray  characteristic,  c,  is  established: 


c  -  n, 


cos 


a0 


(12) 


3 


1.0005  an 
value  of  c 
the  ground 
tions,  are 


typically  is  constrained 
because  the  limits  of  a, 


to  the  interval  of  1.0002  to 
iQ  are  0  to  70  milliradians,  the 
falls  in  the  interval  0.99775  to  1.0005.  Solutions  to 
range,  r,  as  a  function  of  c,  with  the  above  limita- 
presented  in  figure  10.  These  values  were  computed  by 
using  equation  (9),  which  was  numerically  integrated.  As  shown, 
the  quantity  r  varies  from  350  to  500  km  over  the  full  interval 
of  c.  Since  the  discrete  values  of  r,  indicated  by  the  crosses, 
are  well  behaved,  a  fourth  degree  polynomial  interpolation 
function  was  made  that  describes  r  in  terms  of  c.  To  improve 
accuracy,  the  interpolation  function  was  not  fit  directly  to  c 
but  rather  to  a  quantity  defined  as  x  =  (c  -  1)  X  10°  +  2501  . 
The  sum  of  the  square  of  the  deviations  is  less  than  200  m. 
Coefficients  of  the  polynomial  as  fit  to  x  are  listed  in  figure 
10.  The  net  interpolation  function  provides  a  very  rapid  means 
of  computing  r,  as  the  ray  travels  through  the  exponential 
region,  from  knowledge  of  the  surface  refractivity  and  eisvation 
angle . 


CONSTANT  REGION  RAY  TRACING 

The  refractive  index  for  heights  exceeding  50  km  is  assumed 
to  be  unity.  Since  the  refractivity  of  the  medium  does  not  vary 
with  height,  a  ray  propagating  through  this  region  will  not  be 
refracted  and  will,  in  fact,  describe  a  straight  line.  Therefore 
a  geometrical  solution  to  equation  (9)  is  readily  attainable  and 
is  illustrated  in  figure  11.  From  the  law  of  sines, 
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e  =  TT  -  a 

v  v 


arcsin 


~  sin  (<*v  +  n/2) 


(13) 


where  is  the  vector  from  the  earth's  center  to  a  height  of  50 
km,  av  is  the  angle  between  the  normal  of  pv  and  the  ray,  and  ps 
is  the  vector  from  the  earth's  center  to  the  satellite.  If  the 
same  reasoning  is  followed  as  in  the  previous  section,  the  pene¬ 
tration  angle,  ay,  can  be  found  from  Hq  and  by  the  relation¬ 
ship 


=  arctan 


±Apv/z)2-  c2 


(14) 


For  a  satellite  at  1000  km  altitude,  the  rate  of  change  of  r 
with  respect  to  ps  is  approximately  1.46  km  per  km  for  the  inter¬ 
val  0.99775  <  c  <  1.0005.  The  ground  range,  r,  is  normalized  to 
a  satellite  height  of  1000  km  by  adding  (a  +  1000  -  p$)  x  1.46  to 
the  observed  ground  range.  Normal i za t ion  allows  the  library  data 
to  be  computed  only  once,  saving  a  considerable  amount  of  time. 

A  FORTRAN  program  implementing  the  previous  equations  was 
written  for  a  Data  General  NOVA  BOO  minicomputer  with  a  hardware 
floating  point  processor.  All  computations  are  done  with  15 
digit  accuracy,  and  the  ground  range,  r,  for  a  ray  from  the 
surface  to  a  height  of  1000  km,  is  found  in  approximately  100  ms. 


TRANSMITTER 

Satellite  P76-5  was  selected  as  the  rf  source  for  the 
measurement  tests.  Its  orbit  is  sun- synchronous,  which  allows 
observations  to  be  made  at  nearly  the  same  time  each  day.  The 
satellite  was  designed  and  constructed  by  Stanford  Research 
Institute  to  support  the  Defense  Nuclear  Agency  wideband  satel¬ 
lite  experiment  (ref  10).  It  was  placed  into  a  1000  km  orbit  for 
the  purpose  of  carrying  a  multifrequency  crystal-stabilized 
coherent  rf  transmitter  to  characterize  transionospher ic  propaga¬ 
tion  effects.  Ten  frequencies  are  radiated:  one  at  137.7  MHz, 

seven  in  the  378  to  448  MHz  band,  and  one  each  at  1239  and  2891 
MHz.  The  antennas  are  nadir-directed  and  right  circular  polar¬ 
ized,  Effective  radiated  power  is  22  dBm  at  1239  MHz  and  19  dBm 
at  2891  MHz.  In  addition  to  these  desirable  hardware  character¬ 
istics,  ephemeral  data  could  be  obtained,  on  a  weekly  basis,  from 
the  Naval  Astronautics  Group,  Pt  Mugu,  CA. 


To~.  Early  Results  from  the  DNA  Wideband  Satellite  Experiment- 
Complex  Signal  Scintillation,  EJ  Fremouw,  Radio  Science,  vol  13, 
no  1,  p  167-187,  January-Febr uary  1978. 
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RECEIVER 

The  key  to  the  success  of  the  measurement  program  is  the 
ability  to  detect  the  transmitted  signals  as  the  satellite  is 
rising  or  setting  on  the  horizon.  This  region  is  crucial  since 
the  dominant  refractive  effects  occur  for  elevation  angles  less 
than  20  milliradians.  Also  it  is  the  region  where  the  received 
signals  are  lowest  in  intensity.  Other  design  factors  are  the 
accuracy  in  predicting  the  satellite  position  and  velocity 
vectors  from  available  ephemeral  data  and  an  option  to  acquire 
data  from  satellites  other  than  P76-5.  In  accordance  with  the 
above  points,  the  design  goals  are  as  follows: 

Instantaneous  signal  detection 

Dual  frequency 

-145  dBm  sensitivity 

20  dB  dynamic  range 

10  kHz  acquisition  bandwidth 

100  mrad  antenna  bearnwidth 

Real-time  processing 

Balancing  the  design  criteria  led  to  the  configuration  shown 
in  figure  12.  The  high  gain,  wide  bandwidth  rf  preamplifier  sets 
the  overall  noise  figure  to  4.5  dB.  To  attain  the  required  sen¬ 
sitivity,  a  final  IF  bandwidth  of  20  Hz  was  selected.  A  narrow 
bandwidth  such  as  this  combined  with  the  instantaneous  signal 
detection  requirement  abrogates  the  use  of  a  single  IF  filter.  A 
high  speed  digital  signal  processor  (DSP — Spectral  Dynamics  model 
SD-360)  was  chosen  to  implement  512  contiguous  IF  filters  in  the 
range  from  dc  to  10  kHz.  This  unit  is  the  heart  of  the  receiver, 
since  it  provides  a  dual-channel  processor  (two  independent 
filter  banks),  70  dB  dynamic  range,  10  kHz  acquisition  bandwidth, 
and  real-time  processing. 

Instantaneous  signal  detection  is  accomplished  by  downcon¬ 
verting  the  rf  signals  to  a  nominal  IF  of  5  kHz.  Uncertainties 
in  the  satellite  position  and  frequency  drift  limit  precise 
knowledge  of  the  transmitted  signal  to  approximately  ±2  kHz; 
therefore  the  received  IF  falls  between  3  and  7  kHz.  Since  this 
is  well  within  the  10  kHz  acquisition  bandwidth,  the  received 
signal  is  immediately  detectable.  As  the  satellite  moves  through 
its  orbit,  the  radial  velocity  varies,  causing  a  slowly  changing 
Doppler  shift.  To  keep  the  signal  in  the  acquisition  band,  the 
first  local  oscillators  are  manually  adjusted  to  drive  the  IF  to 
its  nominal  5  kHz.  Both  first  local  oscillators  (Hewlett  Packard 
models  8614  and  8616  signal  generators)  are  phase  locked  frequen¬ 
cies  to  1.995  MHz  offset  from  their  respective  transmit  signals. 
The  second  local  oscillator  (Hewlett  Packard  model  5110B/5100A 


synthesizer)  is  set  to  2.0  MHz  and  drives  both  second  stage 
mixers,  generating  the  5  kHz  IF.  An  analog  tape  records  the  IF 
as  a  backup  against  computer  failure. 

Data  processing  and  archiving  are  performed  by  a  Data 
General  NOVA  3/12  computer.  The  computer  controls  the  DSP  and 
sets  it  into  a  dual  channel  FFT  processing  mode  for  analyzing  a 
10  kHz  band.  Data  gathering  is  initiated  by  an  operator  command 
and  normally  begins  when  a  signature  is  observed  on  a  spectral 
display.  From  this  point,  the  computer  logs  the  time  of  day  and 
examines  the  dual  spectrum  every  50  ms  (DSP  analysis  period).  In 
each  period,  the  amplitude  and  frequency  information  of  the 
received  signals  is  stored  for  further  processing.  The  examin¬ 
ation  continues  until  halted  by  an  operator  command. 

Table  1  lists  the  power  budget  factors  for  the  P76-5  trans¬ 
mitter.  Free  space  path  loss  is  the  attenuation  in  radiated 
power  due  to  geometrical  spreading  and  is  given  by  the  summation 

FSPL  =  32.45  +  20  log  f  +  20  log  d,  (15) 


where  f  is  the  frequency  in  MHz  and  d  is  the  distance  from  the 
source  in  km.  The  values  listed  are  computed  for  a  range  of  4000 
km.  Antenna  gains  of  20  and  32  dB  were  measured  by  using  a 
single  ground  station  circular  parabolic  antenna  2  m  in  diam¬ 
eter.  The  antenna  is  fitted  with  a  broadband  horizontally 
polarized  log  periodic  feed  that  passes  signals  from  800  MHz  to 
10  GHz.  Summing  the  ERP,  FSPL,  gain  and  loss  in  table  1  gives 
received  signal  power  levels  of  -127  dBm  (1239  MHz)  and  -125  dBm 
(2891  MHz)  which  must  be  detected  by  the  receiver.  Measured 

receiver  sensitivities  at  the  two  frequencies  are  -144  and  -143 
dBm,  and  the  corresponding  expected  signal- to-noise  ratios  thus 
are  +17  and  +18  dB.  Constructive  interference  of  the  direct  ray 
and  the  sea-reflected  ray  increases  the  received  power  by  6  dB . 
Therefore,  the  dynamic  range  between  maxima  and  minima  is  limited 
to  +23  dB  at  1239  MHz  and  +24  dB  at  2891  MHz. 


EXPERIMENTAL  MEASUREMENTS 

RADIO  MEASUREMENTS 

Eleven  separate  observations  of  the  1239  and  2891  MHz  emis¬ 
sions  of  satellite  P76-5  were  made  during  July  and  August  1978  . 
These  observations  were  conducted  by  NOSC  personnel  at  building 
T323  ,  which  is  situated  at  32°42l  N  latitude  and  117°15’  W  longi¬ 
tude.  The  receiving  antenna  is  35  m  above  mean  sea  level  (msl) 
and  has  an  unobstructed  view  of  the  oceanic  horizon  for  azimuthal 
angles  of  190°  to  340°.  This  antenna  position  allowed  observa¬ 
tions  of  P76-5  while  it  ascended  and  descended  at  the  horizon  on 
the  same  orbit.  Table  2  illustrates  the  geometrical  and  temporal 
characteristics  of  each  observation.  The  column  labeled  Start 
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Power  Level 


Frequency 

1239  MHz 

2891  MHz 

ERP  (dBm) 

+  22 

+  19 

Free  space  path  loss  (dB) 

-166 

-173 

Antenna  gain 

+  20 

+  32 

Polarization  loss  (dB) 

-3 

-3 

Received  signal  level 
(dBm) 

-127 

-12  3 

Receiver  sensitivity 
(dBm) 

-144 

-143 

SNR  (dB) 

+  17 

+  18 

Table  1.  Power  budget  for  Wideband  Satellite  (P76-5) 


Table  2.  PRISM  /measurement  periods. 


Time  (GMT)  indicates  the  time  at  which  the  observation  began. 
The  termination  of  the  measurement  normally  occurred  approximate¬ 
ly  4  minutes  after  this  time.  Average  height  of  the  satellite 
above  msl  is  computed  from  ephemeral  data  for  each  observation 
and  is  listed  in  the  column  labeled  Nominal  Height  (km) .  Each 
measurement  is  indicated  as  ascending  or  descending  at  the  start 
time,  and  the  azimuthal  angle  from  the  receiver  to  the  horizon 
point  is  shown.  The  eleven  observations  were  made  on  six  separ¬ 
ate  days  at  approximately  2000  GMT  (12  noon  local  time);  and 
except  on  2  August,  both  the  ascending  and  descending  measure¬ 
ments  were  taken  within  20  minutes  of  each  other. 


REFRACTIVE  PROFILE  MEASUREMENTS 


Upper  air  soundings,  known  as  radiosonde  observations  or 
RAOBs,  were  made  by  NOSC  personnel  at  approximately  the  same  time 
as  each  day's  period  of  radio  measurements.  The  radiosondes  were 
launched  at  a  site  only  tens  of  metres  distant  from  the  receiving 
antenna.  To  provide  a  better  definition  of  the  refractivity 
structure , the  humidity  and  temperature  elements  were  electrically 
interchanged.  This  modif ication  allows  the  radiosonde  to  sample 
and  transmit  humidity  information  for  a  longer  period  and  does 
not  affect  the  calibration  or  analysis  of  the  data.  Telemetry 
from  the  radiosonde  was  recorded  on  strip  charts  and  manually 
decoded  to  provide  values  of  pressure,  air  temperature,  and 
humidity.  These  atmospheric  quantities  were  then  processed  by  a 
computer  program  to  obtain  the  modified  refractivity,  in  M-units, 
as  a  function  of  height  in  metres  above  msl.  A  description  of 
the  observed  refractivity  profiles  is  provided  in  table  3.  This 
table  shows  whether  the  duct  is  surface  based  (SBD)  or  elevated 
(ELV)  and  lists  the  geometrical  characteristics.  Strong  surface- 
based  ducts  were  observed  on  20  July  and  16  August;  moderate  to 
weak  elevated  ducts  were  ob served  for  the  remaining  four  periods. 


MEASUREMENT  PERIODS 

A  summary  of  the  inferred  data  is  shown  by  table  4.  The 
summary  includes,  for  each  frequency,  the  inferred  duct  type,  the 
standard  deviation  of  the  range  differences  between  the  observed 
and  "best  fit"  null  locations,  and  a  subjective  measure  of  the 
agreement  between  the  observed  and  inferred  refractivity  pro¬ 
files.  A  rating  of  good  is  given  if  the  inferred  profile  geome¬ 
try  reasonably  approximates  the  observed;  moderate  indicates  that 
the  inferred  duct  type  agrees  with  the  observed  but  the  geomet¬ 
ries  differ  signf lean tly ;  poor  is  assigned  if  neither  the  geomet¬ 
ry  nor  the  type  agrees.  For  the  1239  MHz  data,  the  summary  shows 
that  the  inferred  duct  type  (surface  based,  elevated,  or  no 
ducts)  agrees  with  the  observed  in  seven  of  the  eleven  measure¬ 
ments.  But  the  subjective  agreement  between  the  inferred  and 
observed  is  rated  as  good  in  only  four  cases.  The  inferred  duct 
type  for  the  2891  MHz  data  matches  in  five  of  the  eleven  observa¬ 
tions  although  the  subjective  agreement  of  the  inferred  to  the 
observed  profile  is  moderate  at  best. 


Plots  of  the  received  signal  and  the  inferred  profiles  are 
included  for  reference.  Each  of  the  eleven  measurement  periods 
(all  in  1978)  is  examined,  and  the  observed  refractivity  profile, 
signal  received,  and  inferred  profile  are  discussed. 

20  JULY  MEASUREMENTS 

This  measurement  is  characterized  by  an  intense  surface- 
based  duct,  illustrated  in  figure  13.  Top  height  of  the  duct  is 
359  m  and  the  layer  thickness  is  10  m.  The  M-unit  difference 
between  the  top  and  base  of  the  layer  of  -72.6  M-units  shows  this 
to  be  one  of  the  most  intense  ducts  yet  observed  in  the  San  Diego 
area.  An  interesting  feature  of  this  profile  is  its  obvious 
trilinearity.  Since  the  comparison  library  profiles  are  linear 

segments  (to  9  km),  a  reasonably  good  agreement  is  expected. 

Figure  14  is  a  plot  of  the  1239  MHz  signal  received  as  P76-5 
was  ascending  at  1946  GMT.  The  "best  fit"  pld  curve  is  superim¬ 
posed  and  indicates  a  surface  based  duct  as  the  pld  curve  termin¬ 
ates  at  null  3.  Crosses  (  +  )  on  this  plot  show  the  null  number 
(ordinate)  vs  range  (abscissa)  matrix  of  all  possible  nulls 

searched.  These  crosses  are  shown  for  only  the  even  null  num¬ 
bers,  since  placing  crosses  for  each  null  number  would  clutter 
the  figure.  The  "best  fit"  modified  refractive  index  profile  is 
shown  in  figure  15.  Note  the  spectacular  similarity  between  the 
observed  and  predicted  profiles  even  though  the  top  height  of  the 
predicted  duct  is  109  m  below  the  observed  duct. 

Figures  16  and  17  illustrate  the  2891  MHz  received  signal 

and  the  inferred  profile  for  the  ascension  of  P76-5.  Note  in 
particular  the  uneven  range  spacing  of  the  possible  null  loca¬ 

tions  and  the  ragged  appearance  of  the  signal.  Although  the 
inferred  profile  is  a  surface-based  duct,  the  geometries  are 

considerably  different  from  the  observed  profile. 

Figure  18  illustrates  the  1239  MHz  signal  received  as  P76-5 
was  descending  to  the  horizon,  taken  approximately  16  minutes 
after  the  ascending  measurement.  Note  that  the  possible  null 
locations  are  not  as  deep  as  the  locations  shown  by  figure  14. 
Because  the  computer  did  not  start  acquisition  of  the  real-time 
signal  properly,  the  data  displayed  by  figure  18  is  derived  from 
reprocessing  the  real-time  signal  captured  on  the  analog  tape 

backup  system.  Since  this  tape  tends  to  increase  the  noise 

threshold  of  the  signal,  the  depth  of  the  nulls  is  not  as  great. 
Nevertheless,  the  null  locations  for  ranges  less  than  3350  km  are 
clearly  defined.  The  inferred  profile  is  plotted  as  figure  19. 
This  profile  indicates  a  moderately  strong  surface-based  duct, 

but  its  similarity  to  the  observed  is  not  as  striking  as  in  the 
ascending  case.  Note  also  that  the  standard  deviation  of  the 
range  differences  is  3.13  km  as  compared  to  a  standard  deviation 
of  1.93  km  for  the  earlier  measurement. 

Figures  20  and  21  illustrate  the  2891  MHz  received  signal 
and  the  inferred  profile  for  the  descending  pass  of  P76-5.  Since 
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the  signal  was  processed  from  the  backup  tape  system,  the  possi¬ 
ble  null  locations  are  not  as  deep  as  expected.  The  inferred 
profile  is  a  surface-based  duct  and  appears  to  better  represent 
the  observed  as  compared  to  the  inferred  profile  from  the  lower 
frequency  measurement.  The  best  fit  standard  deviation  is  2.11 
km  for  the  2891  MHz  signal  and  3.13  km  for  the  1239  MHz  signal. 

28  JULY  MEASUREMENTS 

This  measurement  is  characterized  by  a  moderate  elevated 
duct  and  is  shown  in  figure  22.  Top  height  of  the  duct  is  787  m, 
with  a  layer  thickness  of  321.7  m.  Note  that  several  linear 
segments  compose  the  trapping  layer.  The  base  of  the  duct  is 
266.8  m  and  the  duct  thickness  is  520.2  m. 

Figure  23  is  a  plot  of  the  1239  MHz  signal  received  as  P76-5 
was  ascending  at  1951  GMT.  The  possible  null  locations  are  rela¬ 
tively  well  defined,  with  amplitude  depths  of  nearly  20  dB.  The 
inferred  profile  (fig  24)  shows  a  nonducting  situation,  which 
does  not  agree  with  the  observed.  Standard  deviation  is  5.28  km, 
a  magnitude  which  indicates  that  the  inversion  technique  per¬ 

formed  poorly. 

Figures  25  and  26  illustrate  the  2891  MHz  signal  and  the 

inferred  profile  for  the  ascension  of  P76-5.  The  range  separa¬ 
tion  of  the  possible  null  locations  is  uneven  and  the  appearance 
of  the  signal  is  ragged.  An  intense  surface-based  duct  is  infer¬ 
red.  Surprisingly,  the  standard  deviation  is  1.79  km,  which 

indicates  a  close  approximation. 

Figure  27  is  a  plot  of  the  1239  MHz  signal  received  as  P76-5 
was  descending  at  2006  GMT.  Possible  null  locations  are  well  de¬ 
fined.  An  elevated  duct  is  inferred,  as  illustrated  by  figure 
28.  The  duct  thickness  of  410  m  agrees  favorably  with  the 
observed  thickness,  but  the  top  of  the  duct  is  663  m  higher.  A 
standard  deviation  of  1.41  km  indicates  a  good  fit. 

Figures  29  and  30  illustrate  the  2891  MHz  signal  and  the 

inferred  profile  for  the  descending  pass  of  P76-5.  The  possible 
null  locations  are  unevenly  spaced  and  the  signal  appears  ragged. 
An  elevated  duct  is  inferred  although  the  geometry  is  different 
from  the  observed.  A  standard  deviation  of  2.06  km  suggests  a 
moderately  good  fit. 


2  AUGUST  MEASUREMENTS 

This  measurement  is  characterized  by  a  dual  elevated  duct 
and  is  shown  in  figure  31.  The  lowest  duct  extends  from  a  base 
at  422.5  m  to  a  top  at  693.6  m.  A  second  duct  exists  from  a  base 
at  769  m  to  a  top  at  1171.4  m.  The  observed  structure  cannot  be 
inferred,  since  only  single-duct  structures  are  included  in  the 
library. 
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Figure  32  illustrates  the  1239  MHz  signal  received  as  P76-5 
was  descending  to  the  horizon.  The  ascension  of  P76-5  occurred 
at  an  azimuth  angle  that  the  receiver  could  not  observe.  An  ele¬ 
vated  duct  is  inferred  and,  referencing  figure  33,  is  a  reasona¬ 
ble  one-duct  approximation  to  the  observed.  A  standard  deviation 
of  2.64  km  is  found. 

Figures  34  and  35  illustrate  the  2891  MHz  signal  and  the 
inferred  profile  for  the  descending  pass.  The  relative  spacing 
of  the  possible  null  locations  is  more  regular  although  the 
signal  is  ragged  in  appearance.  An  elevated  duct  is  inferred 
that  is  much  weaker  than  the  observed.  A  standard  deviation  of 
1.38  km  indicates  a  good  fit. 

3  AUGUST  MEASUREMENTS 

This  measurement  is  characterized  by  a  relatively  weak  ele¬ 
vated  duct  illustrated  by  figure  36.  The  duct  extends  from  390.6 
to  584.4  m  with  a  AM  of  -13  M-units.  A  trilinear  profile  could 
be  constructed  that  reasonably  fits  the  observed. 

Figure  37  is  a  plot  of  the  1239  MHz  signal  received  as  P76-5 
was  ascending  at  2020  GMT.  The  range  separation  between  possible 
null  locations  is  uneven  and  the  signal  for  ranges  greater  than 
3200  km  is  ragged.  A  moderate  albeit  thick  surface-based  duct  is 
inferred  and  is  shown  in  figure  38.  The  standard  deviation  of 
5.58  km  indicates  a  poor  fit. 

Figures  39  and  40  illustrate  the  2891  MHz  signal  and  the 
inferred  profile  for  the  ascension.  The  signal  is  typical  of  the 
2891  MHz  data,  having  uneven  spacing  between  the  possible  null 
locations  and  appearing  ragged.  A  standard  deviation  of  2.61  km 
indicate^  a  moderate  fit  although  the  inferred  profile  is  an 
intense  surface-based  duct. 

Figure  41  is  a  plot  of  the  1239  MHz  signal  received  as  P76-5 
was  descending  at  2034  GMT.  Possible  null  locations  are  readily 
observable  for  ranges  less  than  3300  km  and  tend  to  lose  defini¬ 
tion  at  greater  ranges.  An  intense  low  level  surface-based  duct 
is  inferred,  as  shown  in  figure  42.  The  standard  deviation  is 
3.57  km. 

Figures  43  and  44  illustrate  the  2891  MHz  signal  and  the 
inferred  profile  for  the  descending  pass.  The  inferred  profile 
is  an  intense  surface-based  duct  with  a  top  at  750  m.  A  moder¬ 
ately  good  fit  is  indicated  by  a  standard  deviation  of  2.64  km. 


14  AUGUST  MEASUREMENTS 

This  measurement  is  characterized  by  a  moderate  elevated 
duct,  as  shown  in  figure  45.  The  duct  exists  between  128.9  and 
63  3.9  m  with  a  AM  of  -33.5  M-units.  The  profile  is  reasonably 
trilinear. 
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Figure  4b  is  a  plot  of  the  1239  MHz  signal  as  P7b-5  was 

ascending  at  2039  GMT.  Possible  null  locations  are  readily 
observable  for  ranges  less  than  3350  Km  and  tend  to  become  indis¬ 
tinct  at  longer  ranges.  A  moderate  surface-based  duct  is  infer¬ 
red,  as  shown  in  figure  47.  Standard  deviation  is  2.39  km. 

Figures  48  and  49  illustrate  the  2891  MHz  signal  and  the 

inferred  profile  for  the  ascending  pass.  A  moderate  surface- 

based  duct  is  inferred  and  the  standard  deviation  is  2.89  km. 

Figure  50  is  a  plot  of  the  1239  MHz  signal  as  P76-5  was 

descending  at  2051  GMT.  The  signal  is  similar  to  the  ascending 
pass.  A  tremendous  elevated  duct  is  inferred  and  is  shown  in 
figure  51.  Top  height  of  the  inferred  duct  is  1750  m,  with  base 

at  279  m  and  a  AM  of  -90  M-units.  A  standard  deviation  of  3.44 

km  is  found. 

Figures  52  and  53  illustrate  the  2891  MHz  signal  and  the 

inferred  profile  for  the  descending  pass.  A  strong  surface-based 
duct  is  inferred  and  the  base  of  the  layer  closely  matches  the 

observed.  Standard  deviation  is  2.13  km. 


16  AUGUST  MEASUREMENTS 

This  measurement  is  characterized  by  a  strong  surface-based 
duct,  as  shown  by  figure  54.  Top  height  of  the  duct  is  at  394.5 
m,  with  a  AM  of  -58.8  M-units.  The  profile  can  be  accurately 
represented  by  a  trilinear  structure. 

Figure  55  is  a  plot  of  the  1239  MHz  signal  as  P76-5  was 
ascending  at  2013  GMT.  The  possible  null  locations  appear  rela¬ 
tively  clean  and  correctly  separated  in  range.  A  surprisingly 
good  approximation  to  the  observed  profile  is  inferred  as  shown 
by  figure  56,  although  the  top  of  the  duct  is  94  m  below  the 
observed.  The  M  difference  is  -59  M-units  and  closely  matches 
the  observed.  A  good  fit  is  indicated  by  a  1.73  km  standard 
deviation . 

Figures  57  and  58  illustrate  the  2891  MHz  signal  and  the 

inferred  profile  for  the  ascending  pass.  The  possible  null  loca¬ 
tions  are  unevenly  spaced  and  the  signal  amplitude  appears  rag¬ 
ged.  A  strong  elevated  duct  is  inferred  with  a  standard  devia¬ 
tion  of  2,53  km. 

Figure  59  is  a  plot  of  the  1239  MHz  signal  as  P76-5  was 

descending  at  2028  GMT.  The  possible  nulls  are  not  as  clean  as 
in  the  ascending  pass.  Nevertheless,  the  inferred  profile 
reasonably  approximates  the  observed,  as  shown  by  figure  60.  A 
standard  deviation  of  2.83  km  is  found. 

Figures  61  and  62  illustrate  the  2891  MHz  signal  and  the 

inferred  profile  for  the  descending  pass.  A  strong  surface-based 
duct  is  inferred  and  a  standard  deviation  of  2.77  km  is  found. 


RESULTS 


This  techniaue  was  appl ied  to  eleven  observations  of  dual¬ 
frequency  transmissions  made  on  a  sa tel  1 i te- to-ground  path.  The 
major  findings  are  as  follows: 

Seven  observations  of  the  1239  MHz  data  correctly  predicted 
the  measured  ducting  environment.  However,  the  qeometry  of 
the  inferred  profile  (top  height  of  the  duct,  base  of  the 
layer,  etc)  can  differ  s i on i f i can tl y  from  the  measured 
profile  qeometry. 

Five  observations  of  the  2P91  MHz  data  correctly  predicted 
the  measured  ducting  environment.  Aaain,  the  qeometries  of 
the  observed  and  inferred  profiles  can  differ  significantly. 

Ray  traces  (based  on  the  measured  re f r ac t i v i t y )  compared  to 
the  observed  radio  data  indicate  rms  elevation  angle  errors 
of  R0S  urad  for  the  123°  MHz  data  and  ^30  urad  for  the  2P91 
MHz  da ta  . 

The  rms  range  error  between  the  computed  null  location 
(based  on  the  measured  ref ract ivi ty)  and  the  observed  null 
location  is  ^.2  km  for  the  1239  MHz  data  and  3.7  km  for  the 
2P  °  1  MHz  data. 

The  rms  elevation  angle  and  range  errors  appear  to  be  caused 
by  small  scale  refractivity  perturbations.  These  errors  are 
uncor rectabl e  and  bound  the  attainable  accuracy  for  any 
shallow  angle  rf  measurements . 

The  first  two  findings  are  derived  directly  from  table  4, 
which  is  a  summary  of  the  data  inferred  from  the  rf  measurements. 
A  detailed  examination  of  each  measurement  period  is  provided  in 
the  prior  section. 

Range  errors  are  estimated  by  computing  the  null  locations 
by  means  of  the  observed  refractive  structure  and  comparinq  these 
locations  to  the  observed  locations.  The  comparison  is  identical 
to  that  described  for  the  inference  technique.  Fiqures  ^3  and  ^ 
illustrate  the  range  errors  for  all  measurement  periods  as  a 
function  of  elevation  angle.  The  rms  errors  are  ^.2  km  for  the 
1  239  MHz  data  and  3.7  km  for  the  2P1  MHz  data.  Since  the  range 
separation  between  higher  elevation  anale  nulls  is  approx ima tel  y 
20  km  (1239  MHz)  and  P.7  km  (2P91  MHz),  the  inference  technique 
performs  better  with  the  lower  frequency  data  than  with  the 
higher  frequency  data.  An  rms  error  of  3.7  km  is  nearly  half  of 
the  range  separation  at  the  higher  frequency.  Therefore  the  null 
locations  are  uncertain  to  plus  or  minus  one  null  number.  Since 
the  rms  range  error  for  the  lower  frequency  is  approx ima t el y  20% 
of  the  range  separation,  the  uncertainty  in  the  null  location  is 
less  than  one  null  number  and  can  be  extracted  with  better 
rel iabi 1 i ty . 


Elevation  angle  errors,  shown  by  figures  65  and  66,  are  com¬ 
puted  from  the  range  error  measurements.  For  each  null  location 
in  the  observed  data,  an  equivalent  elevation  angle  is  interpo¬ 
lated  from  the  predicted  ray  trace  data.  The  elevation  angle 
error  is  the  difference  between  the  predicted  angle  and  the 
interpolated  angle.  The  0.61  mrad  error  for  the  1239  MHz  signal 
is  in  close  agreement  to  a  value  of  0.52  mrad,  found  by  Crane 
(ref  11)  at  an  elevation  angle  of  17  mrad  for  a  1295  MHz  radar. 


CONCLUSION 

A  complex  measurement  program  was  undertaken  to  examine  the 
feasibility  of  passively  inferring  the  tropospheric  refractive 
index  structure  from  observations  of  low  angle  satellite- to- 
ground  rf  transmissions.  The  results  of  this  study  show  moderate 
success  in  profile  inference,  but  the  technique  cannot  provide 
reliable  and  operationally  significant  data. 


RECOMMENDATION 


Although  the 
in  inferring  the 
recommended  that 
significant  data 


results  of  this  study  show  a  moderate  success 
refractive  structure  from  rf  measurements,  it  is 
this  project  be  terminated,  since  operationally 
cannot  be  extracted  reliably. 
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Figure  2.  Received  signal  as  a  function  of  range  for  a  standard 
atmosphere  and  a  surface-based  duct. 
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Figure  4.  Nonducting,  elevated  duct,  and  surface-based  duct 
regimes  for  library  profiles. 
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Figure  12.  Block  diagram  of  PRISM  receiver. 
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Figure  19.  Inferred  profile,  1239  MHz. 
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i  Figure  29.  2891  MHz  signal,  28  July  at  2006  GMT. 
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Figure  30.  Inferred  profile,  2891  MHz 
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1239  MHz  signal,  2  August  at  19 
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Figure  62.  Inferred  profile,  2891  MHz. 
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Figure  65,  Elevation  angle  error,  1239  MHz. 
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Figure  66.  Elevation  angle  error,  2891  MHz. 
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